Universality in short-range Ising spin glasses 
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The role of the distribution of couphng constants on the 
critical exponents of the short-range Ising spin-glass model is 
investigated via real space renormalization group. A saddle- 
point spin glass critical point characterized by a fixed-point 
distribution is found in an appropriated parameter space. The 
critical exponents /3 and u are directly estimated from the 
data of the local Edwards- Anderson order parameters for the 
model defined on a diamond hierarchical lattice of fractal di- 
mension df = 3. Four distinct initial distributions of coupling 
constants (Gaussian, bimodal, uniform and exponential) are 
considered; the results clearly indicate a universal behavior. 
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The nature of the spin-glass (SG) transition on real 
systems remains a subject of controversies. There is 
no definite agreement among theories, numerical sim- 
ulations and experimental results, about the standard 
critical exponents governing the transition, as it occurs 
for pure and disordered spin systems without frustration. 
Different apparatus used to measure the non-linear sus- 
ceptibilities lUJ^, distinct numerical procedures used in 
computational simulations &fl] and the lack of an ex- 
act renormalization-group scheme to treat short-range 
spin-glass models in finite dimensions {d > 1), produce 
a spread in the values of the critical exponents, lead- 
ing some authors to claim the absence of universality 
in the critical behavior of spin glasses |g]. The most 
popular renormalization-group method used to investi- 



gate spin-glass models on finite hypercubic lattices 0| 
has been the so-called Migdal-Kadanoff scheme B , which 
was proven to be exact, for pure systems, on a class of 
diamond hierarchical lattices (DHL) |9|. Such a tech- 
nique has been valuable in the study of short-range spin 
glasses [[7|jlC|-[l4|] and was used recently to investigate the 
nature of the low temperature phase for an Ising SG 
on a DHL [Q. Usually, to treat the spin-glass prob- 
lem within a renormalization-group (RG) approach one 
has to choose arbitrarily an initial distribution of ex- 
change coupling constants. In the standard approach wj, 
the parameter space of the renormalization flow is one- 
dimensional and the corresponding renormalized quan- 
tity is the mean-square-root deviation of Kij = Jij/ksT, 
where {Jij} are the random quenched nearest-neighbour 
coupling constants, fcs is the Boltzmann constant and T 
is the temperature. For different chosen initial distribu- 
tions P{{Kij}), different critical temperatures and sets 
of critical exponents have been found either using the RG 
scheme or performing numerical simulations P,pl]|. 

In the present work, the role of the initial distribution 
on the critical properties of the Ising spin-glass model 
is investigated. We consider the short-range Ising spin 
glass defined on the DHL with fractal dimension df = 3, 
described by the Hamiltonian, 
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where the cr's represent the Ising variables assigned 
to the lattice sites and the sum is restricted to the 
nearest-neighbour sites. The evolution of the distri- 
butions, under the renormalization process, is anal- 



ysed in an appropriated higher-dimensional parame- 
ter space. A two-dimensional projection of this space 
is the plane < tanh?{Jij / ksT) > versus ksT/ < 



Jt > 



1/2 



In this representation, each symmetric dis- 



tribution P{{Jij}) is assigned to a curve defined by 
FpiT) =< tanh^ (Jij / ksT) >p, where < ... >p means 
an average over the distribution uM- This parameter 
space is separated in two regions, each one governed by 
its respective attractor, the spin-glass (1,0) and param- 
agnetic (0, cxd) stable fixed points. The frontier separat- 
ing these regions is defined by the locus of Fp{T^), T^ 
being the critical temperature, above (below) which the 
fiow evolves to the paramagnetic (spin glass) stable fixed 
point. Within this frontier, one finds the saddle-point 
spin- glass critical point, characterizing the "fixed-point" 
distribution P*({Kij}), which should remain invariant 
under renormalization. In Figure 1, we display a realiza- 
tion of the flow diagram in the vicinity of the saddle-point 
spin-glass critical point, considering four initial distinct 
symmetric distributions, namely the Gaussian, bimodal, 
exponential and uniform distributions. 



total number of sites belong to the last generation, dif- 
ferent values of (3 appear, associated with the behavior 
of the system close to these critical points, suggesting a 
breaking of universality pOJ . 

To overcome this apparent breaking of universality, 
herein we estimate the critical exponents by calculating 
the temperature dependent EA order parameter close to 
the saddle-point spin-glass fixed point, taking into ac- 
count the "fixed-point" distribution P*{{Kij}). Since 
the exact analytical form of this distribution is unknown, 
these exponents can not be determined, unless one care- 
fully probes the critical region. For that, we monitored 
numerically the evolution of the coupling constant dis- 
tribution until it reaches the F* (T) curve at the vicinity 
of saddle-point spin-glass critical point. This is done by 
following the flow in the diagram < tanh^{Jij/kBT) > 
versus ksT/ < Jf, >i/2 and counting the number of steps 
until the fiow reaches the F*{T) curve. For the present 
DHL this behaviour occurs after the fourth renormaliza- 
tion step, as can be seen in Figure 1. 
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From the data of this plot we estimate the critical tem- 
peratures T^ corresponding to each distribution, recov- 
ering some previous results M , as well as the critical tem- 
perature T* corresponding to the "fixed-point" distribu- 
tion as shown in Table 1. We notice that the positions 
of the Fp (T) curves in Figure 1 follow a decrease of the 
kurtosis k =< Jfj > / < Jfj >^ of the corresponding dis- 
tributions, as we go from the left to the right. We have 
estimated the kurtosis of the "fixed-point" distribution, 
obtaining k ~ 3.30, i.e. is close to the Gaussian one. 

In a previous work pQ] , we have made a direct estimate 
of the critical exponents associated with the order param- 
eter (/? exponent) and the correlation length (ly exponent) 
for the corresponding above-mentioned initial distribu- 
tions. In that approach the exponents were computed 
from the set of numerical values of the local Edwards- 
Anderson (EA) order parameter, obtained through an 
exact recursive procedure developed by the authors [|2| . 
Within this methodology, the values of the local EA or- 
der parameter, for the sites introduced at each hierarchy, 
are calculated taking into account the random couplings 
generated by the corresponding renormalized distribu- 
tion. Therefore, in the vicinity and below Tf, the values 
of the local E A order parameter of the last hierarchy were 
calculated with the coupling constants given by the ini- 
tial distribution. Since, for the present DHL, | of the 
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FIG. 1. The evolution of the probability distributions, un- 
der the renormalization process, in a two-dimensional space. 

From this point on, wc apply the method developed in 
]l2[ , calculating the local magnetization of each lattice 
site and the EA order-parameter, varying the tempera- 
ture and taking an average over many initial distribu- 
tion realizations (samples) . The critical exponent associ- 
ated with the EA order parameter per spin is defined by 
QEA ^ [(r* - T)/T*f where T* is the spin-glass critical 
temperature corresponding to the saddle-point spin-glass 
critical point. Moreover, the correlation-length critical 
exponent can be also estimated by scaling the order pa- 



rameter at T = T*, i.e 



EA ~ Mf^/"" where M = 6^ + 1, 
6 being the scaHng factor (herein we restrict ourselves to 
b — 2) and N the number of lattice hierarchies. We con- 
sider the four above-mentioned distributions of coupling 
constants as the initial ones, evaluating the correspond- 
ing /? critical exponents in each case, for lattices up to 
fifteen generations and 200 samples. To obtain the ratio 
P/v we consider N varying from 8 to 15 and the number 
of samples varying from 1000 (for smaller systems) to 300 
(larger systems). 

In Figure 2, we exhibit the plot of the EA order pa- 
rameter per spin as a function of the temperature for 
each case. Note that all plots fall onto the same curve 
within the errors bars. One notices that larger error bars 
are found as one gets away from the critical temperature 
T* , contrary to what happens in conventional numeri- 
cal simulations. This is due to the fact that our sample 
averages are taken over different initial conditions (i.e., 
distinct pools of numbers associated to a given initial 
distribution M). 



generated by the considered algorithm avoid the reach 
of the saturated limit. The corresponding points in the 
flow diagram, indicating where the order parameter was 
calculated, all fall onto the F*{T) curve characterizing 
the spin- glass "fixed-point" distribution, as shown in the 
inset of Figure 2. 

In Figure 3, we exhibit the log-log plot of the data in 
Figure 2. The results from distinct initial distributions 
all coincide within the error bars, indicating a universal 
behaviour for the f3 exponent. The scaling region used for 
the computation of the (3 exponents was ln[{T* — T)/T*] 
in the range from —4.0 to —2.0; small variations in the 
width of this scaling region (e.g., from —4.0 to —1.5) 
did not change significantly our estimates (/3 remainig 
inside the error bars). A similar universal behavior is 
observed in Figure 4, where we exhibit the logarithm of 
Qea as function of the hierarchical level N, at the critical 
temperature T*, yielding the same P/v ratio within the 
error bars. 
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FIG. 2. The EA order parameter per spin as a function 
of the temperature (in units fc_g = 1), in the neighbourhood 
of T* . The inset shows points in the flow diagram were the 
values of qea where calculated. 



FIG. 3. The log-log plot of the EA order parameter per 
spin in terms of the distance from the critical temperature 
Tc . One clearly sees a universal behaviour within the error 
bars. The values of the (3 exponents were computed from the 
points inside the scaling region shown. 



Close to the critical point the distribution parameters 
(e.g., its mean-square deviation) vary slowly and so are 
not very sensitive to such initial conditions, whereas away 
from T*, different sets of numbers may lead to RG tra- 
jectories far apart in the parameters space. For temper- 
atures very close to the critical point, finite-size effects 
occur, preventing the vanishing of the order parameter 
at T*, while for low temperatures, numerical difficulties 
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FIG. 4. The logarithm of the EA order parameter per 
spin as a function of the hierarchical level, at the universal 
spin-glass critical temperature T* . 

In Table 1, we exhibit the corresponding values of /?, 
(3/v and the estimated values of v, corresponding to each 
initial distribution. We call the reader's for the fact that 
these "universal" values for the critical exponents are 
very close to the ones obtained by considering the Gaus- 
sian distribution as the initial one within the approach 
of reference [0. This is an expected result since the 
" fixed-point" distribution is very close to the Gaussian 
one |l^,|ri|]. This can also be seen in Figure 1 by noting 
that the point corresponding to the Gaussian distribution 
at T^ ~ 0.88 is indeed very close to the curve F*{T). 

TABLE L The estimated values of the critical tempera- 
ture T^ for each initial distribution of coupling constants, the 
universal spin-glass critical temperature T* , and the critical 
exponents (5 and f. The columns are disposed according to 
decreasing kurtosis of the corresponding initial distribution. 





Exponential 


Gaussian 


Uniform 


Bimodal 


rr 


0.723 ±0.001 


0.881 ±0.001 


0.998 ±0.001 


1.132 ±0.001 


t; 


0.854 ±0.002 


0.856 ±0.002 


0.856 ±0.001 


0.854 ±0.002 


ii 


0.63 ±0.03 


0.63 ±0.03 


0.63 ±0.03 


0.63 ±0.03 


13/v 


0.36 ±0.01 


0.36 ±0.01 


0.36 ±0.01 


0.35 ±0.01 


V 


1.8 ±0.1 


1.8 ±0.1 


1.8 ±0.1 


1.8 ±0.1 


K 


6.0 


3.0 


1.8 


1.0 



It should be stressed, as discussed in ll^l, that the 
agreement of the present results with those obtained from 
other methods, developed for the cubic lattice (e.g., ex- 
tensive Monte Carlo simulations [Hj|], high-temperature 
series expansions pq]) is very impressive, taking into ac- 
count that the former are, in principle, valid for a, df — 3 
diamond hierarchical lattice. It should be mentioned 
that, as far as the hierarchical lattice is concerned, the 
only approximations involved in the present approach are 
the finite sizes, and the number of samples investigated. 
A further analysis of the good performance of the short- 
range Ising SG model defined on a DHL, with scaling 
factor two and fractal dimension three, was presented in 
ref. in. 

In conclusion, we have shown that, as far as the 
Migdal-Kadanoff renormalization group scheme is con- 
cerned, the critical behavior of the short-range Ising SG 
model is governed by a unique saddle-point spin-glass 
fixed point, characterized by a fixed-point distribution of 
coupling constants. Furthermore, its critical exponents 
are universal and independent of the particular initial dis- 
tribution of coupling constants employed. Whether this 
scenario will persist for more sophisticated approaches 
remains an open question, now under study. 
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